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We establish a direct connection of quantum Markovianity of an open quantum system to its 
classical counterpart by generalizing the criterion based on the information flow. Here, the flow 
is characterized by the time evolution of Helstrom matrices, given by the weighted difference of 
statistical operators, under the action of the quantum dynamical evolution. It turns out that the 
introduced criterion is equivalent to P-divisibility of a quantum process, namely divisibility in terms 
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Moreover, it is shown that similar mathematical representations as those found for the original 
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backflow and establish a local and universal representation of the measure. We illustrate some 
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I. INTRODUCTION 


Even though experimental techniques are developing 
rapidly, perfect isolation of fragile quantum systems 
from noisy environments cannot be warranted in general, 
which necessitates effective descriptions for nonunitary 
dynamics [T]. A well-known treatment of the dynam¬ 
ics of such systems is provided by a quantum dynami¬ 
cal semigroup represented by a generator of the Gorini- 
Kossakowski-Sudarshan-Lindblad form nig. 

However, this description uses several rather drastic 
approximations which do not apply to open quantum 
systems in general. The mismatch is typically asso¬ 
ciated to the neglect of memory effects. Classically, 
there exists a well-established mathematical theory deal¬ 
ing with stochastic dynamics featuring memory effects, 
based on the theory of stochastic processes, yet the def¬ 
inition of classical Markovian stochastic processes can¬ 
not be straightforwardly carried over to the quantum 
regime. This led to an intense debate along with differ¬ 
ent proposal for the characterization and quantification 
of memory effects in the dynamics of open quantum sys¬ 
tems. Among others, approaches to non-Markovianity 
were based on deviation from semigroup dynamics [1], 
on divisibility in terms of completely positive maps [S], 
on dynamics of entanglement [5] and correlations [B] , and 
on the Fisher information [3, see [S] for a review. 

In this work we focus on a straightforward extension 
of the measure of non-Markovianity introduced in Ref. 
[iiin] and reviewed in m, which characterizes memory 
effects by an exchange of information between the open 
quantum system and the environment. This interpreta¬ 
tion applies also to the generalized measure that uses Hel¬ 
strom matrices known from quantum estimation theory 


|12] . A similar approach has already been introduced in 
Ref. [13] focusing, at variance with the present proposal, 
on divisibility in terms of completely positive maps (CP- 
divisibility). However, CP-divisibility yields only a suffi¬ 
cient condition for quantum Markovianity with respect to 
the novel characterization proposed in this article, which 
turns out to be equivalent to divisibility in terms of posi¬ 
tive maps (P-divisibility), and has in addition a clear-cut 
connection to classical Markov stochastic processes. The 
generalized criterion introduced in Section ^ thus com¬ 
bines a physical interpretation of non-Markovianity in 
terms of an information flow and a relation to the well- 
known classical definition. 


Moreover, the associated measure has similar mathe¬ 
matical features and representations as those found for 
the original definition miiB] which simplify its analyt¬ 
ical, numerical and experimental determination drasti¬ 
cally. Indeed we prove in Section HI A that, first, optimal 
initial states for non-Markovian dynamics, experiencing 
a maximal backflow of information, must be orthogonal 
and, second, the measure admits a local representation 
showing locality and universality of quantum memory ef¬ 
fects. 


An illustration of the essential features of the gener¬ 
alized characterization of non-Markovianity in compari¬ 
son with the original definition and approaches based on 
CP-divisibility is provided in Section |lHB[ The examples 
show the sensitivity of the novel approach to memory ef¬ 
fects of dynamics to which the original definition was un¬ 
susceptible and illustrate the existing difference between 
different types of divisibility of a dynamical process. Fi¬ 
nally, we summarize our results in Section |TVj 
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II. STATE DISCRIMINATION AND 
P-DIVISIBILITY 

Henceforth, T-L refers to the Hilbert space of the open 
quantum system and the corresponding set of physical 
states, i.e., the set of positive trace class operators with 
unit trace, is denoted by S{'H). Moreover, we assume 
that the dynamics of the open quantum system is deter¬ 
mined by a one-parameter family d) = {<i)t|0<t<T} 
of completely positive, trace preserving linear maps <i>t, 
where <i)o refers to the identity [T]. 

A. Ensemble discrimination 

The previously introduced definition of quantum non- 
Markovianity laiin], which was based on the concept 
of an information flow between system and environment, 
relied on the trace distance between two quantum states 
Pi and p 2 defined by m 

'Dipi, P 2 ) = 1\\P1 - P2\\l ■ (1) 

Here ||H||i = Tr|A|, where the modulus of an operator 
is given by |A| = -s/aTA, refers to the so-called trace 
norm on the set of trace class operators m- A natural 
and interesting generalization of this quantity is obtained 
when one also allows for biased probability distributions. 
That is, one considers 

IIpiPi-P 2 P 2 II 1 = Tr|pipi-P 2 P 2 I , (2) 

where {pi} refers to an arbitrary binary probability dis¬ 
tribution, i.e. pi _2 > 0 and pi +P 2 = 1- Clearly, choosing 
an unbiased distribution, i.e. pi ^2 = 1/2, one obtains 
the original expression Q. This generalization was first 
proposed in this context by Chruscihski et al. ca. 

The Hermitian operator A = PiPi —P 2 P 2 is also known 
as Helstrom matrix and we find, applying the 

triangle inequality to Eq. Q, 

bi-P 2 I < llA|li < Pi-kp2 = 1 • (3) 

Clearly, the lower bound is obtained for pi = p 2 whereas 
the upper bound is attained if and only if pi T p 2 mean¬ 
ing that the two states have orthogonal support (which is 
defined as the subspace spanned by the eigenvectors with 
nonzero eigenvalues). This is easily shown by means of 
the characterization of the trace norm 

||A||i = 2maxTr{nA}-I-P 2 -Pi , (4) 

where the maximum is taken over all projection operators 
H. Here, Eq. ([^ is derived employing the Jordan-Hahn 
decomposition of the Hermitian Helstrom matrix A in 
terms of two positive and orthogonal operators, i.e. A = 
S — Q where S,Q >0 and S ± Q. 

The interpretation of the trace distance of two states 
Pi ,2 S S{'H) as a measure for their distinguishability di¬ 
rectly carries over to the trace norm of a Helstrom ma¬ 
trix A. Consider a one-shot, two state discrimination 


problem where Alice prepares one out of two quantum 
states pi ^2 with corresponding probability pi, 2 , he., we 
have Pi -l- P 2 = 1- Hence, we allow for the general situ¬ 
ation of a biased preparation by Alice. She finally sends 
the prepared state to Bob who performs a single (strong) 
measurement to infer which state he had received m- 
To guess the state from the measurement with possible 
outcomes ft, Bob defines two sets of possible results, R 
and H\i?, and assigns the state to be pi if the measure¬ 
ment outcome is in R and p 2 if a value in f2\i? is ob¬ 
tained. This strategy results in an effective two-valued 
positive operator valued measure (POVM) {Tr, 1 — Tr} 
where Tr refers to the collection of effects correspond¬ 
ing to outcomes in R. The probability for correct state 
discrimination via this strategy is then given by 

^success = piTrjTRpi} -I- P2Tr{(l - rR)p2} 

= P2 -I- Ti'IATr} (5) 

which is maximal if Tr is the projection n{A>o} on 
the subspace spanned by eigenvectors corresponding to 
positive eigenvalues of A. Employing that Tr|A| = 
Tr{A(n{jf>o} — n{jc<;o})} for any Hermitian operator 
X one shows that the maximal success probability for 
correct discrimination obeys m 

process = o^ax^T.uccess = ^ (^ + H^Hl) ' (6) 

Hence, the trace norm of the Helstrom matrix A = 
PiPi —P 2 P 2 is the bias in favor of the correct state identi¬ 
fication of the state prepared by Alice so that it may be 
interpreted as a measure for the distinguishability of the 
two states pi and p2 with associated weights pi and P2- 

Due to a result by Kossakowski [HHO] the trace norm 
can also be used to witness positivity of a trace preserving 
map A. That is, a trace preserving map A is positive if 
and only if it defines a contraction for any Hermitian 
operator X with respect to the trace norm, i.e., 

|!AA||i < ||A||i , for all A = At. (7) 

We thus obtain for two states evolving according to a 
dynamical map <i>( 

\\pi<^t{pi) -P2<^t{p2)\\i < WpiPi -P 2 P 2 II , (8) 

for all t where we also used linearity of the map. 

Adopting the previous characterization for quantum 
non-Markovianity laiin] one then directly derives a gen¬ 
eralized criterion which still relies on the concept of an 
information flow: 

Definition 1 . A quantum process d* is defined to be 
Markovian if Wpi^tipi) — P2^t{p2)\\i is a monotonically 
decreasing function of t > 0 for all sets {pi,pi} with 
Pi > 0 , Pi + P2 = ^ and Pi G S{'H). 

We stress that the authors of Ref. m introduced a 
different definition for Markovianity, which used the Hel¬ 
strom matrices of a dilated system making it equivalent 
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to CP-divisibility of the dynamical process. However, as 
we shall see in Section II B[ Markovianity as defined here 
can be related to the concept of P-divisibility and pro¬ 
vides a clear-cut connection to classical Markov stochas¬ 
tic processes. 

We conclude this section by showing that the previ¬ 
ously developed explanation to substantiate the inter¬ 
pretation of quantum memory as an information back- 
flow from the environment to the system m also applies 
here. To this end, we define the quantities 


lint(i) = \\PlPs\t) -p2Ps\t)h 


(9) 


and 


^'extW = Ibipy^(t) -p2Psl(t)||l , (10) 

where Pse(^) ~ UtP^g ®PeUI refers to the state of system 
and environment at time t subject to a unitary evolution 
Ut- We thus have Ps\t) = TrEPg\;{t) and, similarly, the 
state of the environment at time t is given by = 

As discussed in m iZiint(t) describes the distinguisha- 
bility of the open system at time t, while Xext(i) can be 
interpreted as the information on the total system which 
is not accessible by measurements on the open system 
only. Due to contractivity of the trace norm under posi¬ 
tive and trace preserving maps (cf. Eq. 0 ) we conclude 
that both quantities are positive as the partial trace is 
(completely) positive and preserves the trace. Moreover, 
existence of the quantum dynamical map implies that the 
initial state is factorized, therefore Xext(O) = 0, so that 
unitary invariance of the trace norm implies 

Ii„t(t) -k lext(t) = 2:int(0) = const . (11) 

This relation clearly expresses the idea of exchange of 
information between the open system and the environ¬ 
ment as an increase of 2'ext(f) necessarily forces Xint(f) 
to decrease. One may also derive upper bounds for the 
external information which, however, do not allow for a 
clear interpretation in terms of formation of correlations 
between system and environment and changes in the en¬ 
vironmental states as proven for the trace distance in Ref. 

m- 


B. Distinguishability and divisibility 

In this section we introduce the concept of divisibility 
of a dynamical map which has been at the basis of several 
other proposals for quantum non-Markovian dynamics [5J 
Eg and elucidate its connection to the distinguishability 
measure defined before. 

We start by revisiting the notion of n-positivity. A map 
A on "H is said to be n-positive if and only if the linear 
tensor extension of the map given by (A 0 1„)(A ^ B) = 
A(A) (g) B for linear operators A on B and B on C” is 


positive, that is, (A (g 1„)C > 0 for all positive linear 
operators C on H 0 C”. 

Clearly, 1-positivity is equivalent to the map A be¬ 
ing positive and the concept of n-positivity is hierarchi¬ 
cal, which means that n-positivity of A implies that the 
map is also fc-positive for all 1 < fc < n m- The con¬ 
verse, however, is not true in general mug- Finally, one 
speaks about complete positivity if a map is n-positive 
for all n € N. For a finite-dimensional system with 
dimH = N this property is equivalent to Wpositivity 

[ig. 

If the maps <i>t comprising a dynamical process are 
invertible for all t > 0 with inverse then one can 

define a two-parameter family of maps given by 

, ( 12 ) 


for all t > s > 0 such that = ‘I’t and <!>* = 

We remark that although the existence of a left-inverse 
requiring injectivity of the maps would be sufficient, 
we assume to be the left- and right-inverse. The 
notion of divisibility deals with the properties concerning 
n-positivity of these maps, i.e.. 

Definition 2 . A dynamical process <i> is called P-divisible 
(CP-divisible) if ^t,s is a positive (completely positive) 
map for all t > s > 0 . 

Even though is completely positive for all times by 
definition, this does not imply (complete) positivity of 
in general as the inverse of a CP-map need not be 
positive showing that the above definition is nontrivial. 
However, the concept of divisibility relies on the existence 
of the two-parameter family which is limited to processes 
for which exists for all times. This property though 
cannot be taken for granted as e.g. the damped Jaynes- 
Cummings model on resonance EJEHlin] or examples on 
quantum semi-Markov processes |27j show, thus making 
the concept of divisibility sometimes ill-defined. 

Typically, however, the inverse maps exist apart from 
isolated points in time so that one may describe the dy¬ 
namical process by a time-local master equation for the 
open system during intermediate intervals. Assuming a 
sufficiently smooth time dependence, the generator obeys 
ICt = has the general structure 

ICtPsit) = - i [Hs{f),ps(t)] + '^-fj{t)(Aj{t)ps{t)A'){t) 

3 

-^{A]{t)A,{t),ps{t)}) , (13) 


which is similar to the well-known Lindblad form nig 
apart from the fact that the system Hamiltonian Hs{t), 
the rates jj{t) and the Lindblad operators Aj{t) may de¬ 
pend on time as the process might not represent a semi¬ 
group. 

The maps of the two-parameter family are then given 

by 


$ 


t,s 


T<_ exp 



drlCr 


(14) 
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where T<_ refers to chronological time-ordering. We thus 
have K-s = |t=s so that we conclude: 4) is (C)P- 

divisible if and only if K-t is the generator of a (com¬ 
pletely) positive semigroup for all fixed < > 0. 

The famous Gorini-Kossakowski-Sudarshan-Lindblad 
theorem HE] and a result on generators of positive semi¬ 
groups by Kossakowski m then allow for a characteri¬ 
zation of CP- and P-divisibility in terms of the generator 

m- 


Theorem 1. 


The dynamics generated hy K-t (13) 


1. is CP-divisible if and only if"fj{f) > 0 holds for all 
j and t > 0. 


2. is P-divisible if and only if for all n ^ m 


^7j(t)|(m|Aj(t)|n)p > 0 , (15) 

3 


We may alternatively characterize P-divisibility employ¬ 
ing the contraction property 0 and link it to quantum 
Markovian behavior defined above (see definition . To 
this end, we first prove the following lemma. 

Lemma 1. For any Hermitian operator X 0, there 
exists a real number A > 0 and a Helstrom matrix A 
such that X = AA. 

Proof. Let X = X^ 0 be given. If A > 0, then 
Pi = {TrX)~^X defines a state so that for A = TrA 
the Helstrom matrix characterized by pi = 1, P 2 = 0 
with arbitrary p 2 proves the claim. Similarly for A < 0. 

Hence, suppose that A is indefinite. Employing the 
Jordan-Hahn decomposition we thus find nonzero oper¬ 
ators Fi .2 > 0 where Yi T Y 2 such that A = Yi — Y 2 
and, therefore, Tr|A| = TrYi -|- TrY 2 > 0. Clearly, the 
operators pi = {TrYi)~^Yi define states and we have 


holds for any orthonormal basis {|n)} ofTL and all 
t > 0. 

Proof. The first statement is precisely the Corini- 
Kossakowski-Sudarshan-Lindblad theorem HE] and the 
second can be derived from Kossakowski’s result on gen¬ 
erators of positive semigroups m which states: A dy¬ 
namics generated by C is P-divisible if and only if for 
any set of projections H = {Ilm}Tn£i defining a resolu¬ 
tion of the identity, i.e. = Iw, the following 

relations 


0 , 

TO G J , 

(16) 

^mn(n) ^ 0 , 

m n G I , 

(17) 

^ ^ *3-77177,(11) = 0 , 

n G I . 

(18) 


mGl 


for amn(n) = Tr{nm(£n„)} are satisfied. 

It is easily seen that condition (181 refers to the preser¬ 


vation of trace which is always met for the generator ICt 
(131 by its very structure. Rearranging terms of Eq. (18) 


one then finds 


X = X{pipi - P 2 P 2 ) , (21) 

for A = Ti'lAj and pi = A“^TrYi. These quantities are in¬ 
deed positive and sum to one, thus, representing a prob¬ 
ability distribution which concludes the proof. □ 

Since dynamical maps 4>t as well as the trace norm are 
homogeneous (with respect to positive numbers) it thus 
suffices to apply the characterization of positivity 0 to 
Helstrom matrices which yields: 

Theorem 2. If the dynamical maps defining a process 
4) are bijective, then 4) is Markovian if and only if it is 
P-divisible. 

Proof. We first note that pi^t{pi) — P 2 ^t{p 2 ) = 4>t(A) 
holds for any probability distribution {p^} and pair of 
states Pi. Hence, it suffices to consider the time evolution 
of Helstrom matrices when studying quantum Marko- 
vianity. Suppose 4) is P-divisible. It follows that 

|l4>*(A)||i = ||4>*,,(4>,(A))|U<||4>,(A)|U, (22) 


a„n(n) = - ^ amn(n) , (19) 

m^n 

for all n G I. The constraints thus reduce to the single re¬ 
lation 0 which can additionally be restricted to sets of 
rank-one projections H by virtue of linearity. Evaluating 
«mTt(n) for m ^ n and rank-one projections associated 
to an orthonormal basis {|n)} of TL one finds 


amn(n) = Tr{|m)(m|/Ct(|n)(n|)} 

= ’ ( 20 ) 
3 


which is condition (151. 


□ 


Obviously, the conditions for P- and CP-divisibility co¬ 
incide for a master equation with a single decay channel. 


for alH > s > 0 and Helstrom matrices A due to positiv¬ 
ity of 4)( s. Hence, ||4)((A)||i is a monotonically decreas¬ 
ing function of time for any A showing that the process 
4) is Markovian. 

For the converse, we first note that the inverse map 
4))"^ exists for all t as the dynamical maps are bijective 
on the set of Hermitian operators by assumption. Thus, 
the maps 4)t_s exists for all t > s > 0. Now, let the 
process be Markovian, i.e., ||4)((A)||i < ||4)s(A)||i for all 
t > s > 0 and Helstrom matrices A. We may rewrite this 
as 


||‘i>M(‘J>.(A))lli<ll‘J>.(A)lli ■ (23) 

from which positivity of 4)t ^ follows according to Eq. Q 
employing lemma and the fact that 4)t is bijective for 
all t. Hence, 4) is P-divisible. □ 
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C. Connection to classical Markovian stochastic 
processes 


The definition of Markovianity stated above is thus 
equivalent to P-divisibility of a dynamical process d* if 
this notion is well-defined at all. However, the measure 


for non-Markovianity (371 may be evaluated for any dy¬ 


namical process showing its great benefit. Similarly, one 
could obtain equivalence of CP-divisibility with Marko¬ 
vian behavior if the dilated process is considered 

as shown in Ref. m- However, P-divisible quantum 
processes offer the relevant feature of a distinct connec¬ 
tion to classical Markovian stochastic processes. To show 
this, we employ the characterization of P-divisible pro¬ 
cesses given in theorem]^ Let p{t) be the solution of the 
master equation 


dt 


p{t) = ICtp{t) , 


(24) 


with initial state p(0) where ICt is given by Eq. (13). The 
time-evolved state admits an instantaneous spectral de¬ 
composition. 


p{t) = 


(25) 


so that {\4>m{t))} defines an orthonormal basis on H and 
{Pmit)} represents a classical probability distribution for 
all t > 0. By virtue of the orthonormality the eigenval¬ 
ues Pm{t) = {4>m{t)\p{t)\(j)m{t)) obey the following closed 
differential equation: 


dt 

where 


Pmit) ^Wmnit)Pnit) - VP„m(t)pm(t) , (26) 


Wmnjt) — y^7j(t)|(((>m(^)|7lj(t)|0n(^))| • (27) 


Obviously, the term in Eq. (261 with m = n drops out. 


Given any solution of a quantum master equation, one 
thus obtains a classical jump process. We emphasize that 
the rates (t) depend in general on the initial eigenba- 
sis {|</i’m(0))} as well as the initial probability distribution 

{Pm(0)}. 

According to the theory of classical stochastic pro¬ 
cesses, any hierarchy of n-point probability distributions 
PniUnjtn', • ■ • i ?/ij ^i) with discrete sample space O satis¬ 
fying the consistency relations 


PniUrij ■ 5 G) ^ 0 j 
PiiVi^h) = 1 , 




^ ^ Pnijjn') ■ - ‘i Vm: ^m‘i ■ • ■ : ^1: ^l) 
— Pn—l{yni - ■ ■ ! ^l) ; 


(28) 

(29) 


(30) 


for all 1 < m < n and 0 < G < ^2 < ■ • • < in, determines 
a stochastic process y{t) with values in ft [2HH2n]- Here, 
P-niVn, tn] ■ ■ ■', ?/i, G) gives the probability to observe the 
values Pi at times ti for i = 1,.. . ,n for the stochastic 
process y{t). For classical Markov processes obeying the 
relation 


Pl\n—liyn7 tnlVn — l: tn—1 7 • ■ ■ Pi 7 ti) 

~ Pl\liyn7tn\yn—l7tn—l) 7 (^1) 

for any conditional probability distribution the 

full hierarchy is completely determined by the 1-point 
probability distribution Pi and the so-called transition 
probability Piq. However, these two nonnegative func¬ 
tions cannot be chosen arbitrarily but must satisfy 


Piiy27t2) = Y -Pi|i(2^2,i2|yi,ti)-Pi(2/i,G) , (32) 

yiGSl 

-Pi|i(2/3,i3|yi,G) (33) 

= Y -Pi|i(y3,i3|y2,t2)-Pi|i(y2,t2|yi,P) , 


where Eq. (331 is called the Chapman-Kolmogorov equa¬ 
tion. These equations thus characterize uniquely a 
Markov process. 

One may equivalently write this relation for the tran¬ 
sition probability in its differential form 


^Pi|i(y,t|x,s) (34) 


where the transition probability per unit time Wyzit) is 
nonnegative and represents the probability for a transi¬ 
tion to y given the classical state was z at time t. It is 
clear that the (differential) Chapman-Kolmogorov equa¬ 
tion characterizes the transition probability Piq but a 
similar equation, the so-called Pauli master equation, can 
be derived for the 1-point probability distribution, that 
is. 


(35) 

= E [^yzit)Piiz,t) - Wzyit)Piiy,t) . 


One thus concludes that Eq. (261 can be interpreted 
as Pauli master equation for the 1-point probability 
distribution of a classical Markov process with H = 
{1,..., dimP} if and only if 


W^nit) = Y'tjit)\{(l>mit)\Ajit)\(l)nit))\‘^ > 0 , (36) 

3 

for all t > 0 and m ^ n. Hence, P-divisibility of the 
quantum process is a sufficient condition to warrant pos¬ 
itivity of the rates Wmnit) (see Eq. (El)- This shows 
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that quantum non-Markovianity defined with respect to 
the generalized trace distance based measure does not 
only provide an interpretation in terms of an informa¬ 
tion backflow but also allows for a connection to clas¬ 
sical Markovian stochastic processes. Namely, to each 
P-divisible quantum process, given as the solution of a 
master equation of the form (131 with Lindblad opera¬ 


tors and rates satisfying Eq. (151 for an initial state /9(0), 


one associates a classical Markovian process obtained as 


solution of the classical master equation (261 with transi¬ 


tion rates given by Eq. (27) and initial condition specified 
by the eigenvalues of p(0p 

Finally, we note that P-divisibility of the quantum pro¬ 
cess is indeed equivalent to the positivity of the rates 
Wjnn{t) of the classical Pauli master equation if the quan¬ 
tum master equation has the property that the eigenbases 
{|(^ri(t))} of p{t) run over all orthonormal bases when 
varying the initial state /o(0). This is satisfied if the max¬ 
imally mixed state is in the image of the dynamical map 
<i)t which holds true for sufficiently small times due to 
continuity of the process and as we have <i)o = 1. For 
two-level systems it is shown in Appendix that this 
constraint is not only sufficient but also necessary. If any 
orthonormal basis is encountered all classical processes 
derived from the quantum master equation are classi¬ 
cally Markovian if and only if the quantum dynamics is 
P-di visible. 


III. GENERALIZATION OF THE TRACE 
DISTANCE BASED NON-MARKOVIANITY 
MEASURE 


this feature was only sufficient in the original definition 
(see Section IIIB). 


To begin with, we concentrate on the maximiza¬ 
tion procedure contained in the quantifier for non- 
Markovianity (37) which, contrary to the original defi¬ 


nition [3], now even requires to sample also over binary 
probability distributions. Fortunately, a similar charac¬ 
terization of pairs of states maximizing Eq. (371 can be 


proven m and yet also a local representation is admitted 
m simplifying the sampling significantly. 

We finally illustrate the novel definition of non- 
Markovianity by means of examples which show the dif¬ 
ference to the original characterization and other ap¬ 
proaches to non-Markovianity. 


A. Expressions of the generalized 
non-Markovianity measure 

A set {pi,pi, 2 }, where pi ^2 S S{'H) and {pi} defines 
a binary probability distribution, is said to be optimal if 
the maximum in Eq. (371 is attained for it. Note that the 


quantum states of an optimal set are necessarily nonequal 
as the dynamics of the norm of a nonindefinite Helstrom 
matrix under any dynamical map is trivial due to trace 
preservation and positivity of the map. 

Theorem 3. The states of an optimal set must be or¬ 
thogonal, i.e. 


= max / 

{pi,pi±P2} Ja>0 


dt a{t,pi,pi) . (39) 


Having stated the generalized definition for quantum 
non-Markovianity which allows for a connection to clas¬ 
sical Markov processes, we now consider the expression 
of the corresponding measure and address its mathemat¬ 
ical and physical features. According to the definition 
of a quantum Markovian process given in Sec. IIA it is 
natural to consider the following measure for quantum 
non-Markovianity, quantifying the degree of memory ef¬ 
fects with respect to this generalized definition 

A/'(d') = max / dta{t,pi,pf) (37) 

{Pi,Pi} Jcr>0 


with 


a{t,pi,pi) = ^lbi^'t(pi) -P2^'t(p2)||i , (38) 

where the integration runs over all intervals where the 
distinguishability \\pi^t{pi) — P 2 ^t{p 2 )\\i increases. A 
process is said to be non-Markovian if A/’(d)) > 0. As 
discussed in Sec. Ill Al this measure still admits the inter¬ 
pretation as a quantifier for the information flow from 
the environment back to the open system. In addition, 
as discussed in Sec. |IIB| Markovianity will now be equiv¬ 
alent to P-divisibility of a quantum process. Note that 


Proof. Suppose {pi, pip} be optimal with pi f- p 2 . Along 
the lines of the proof of lemma we obtain a probability 
distribution {qi} and two orthogonal states gi ^2 such that 

PiPi -P 2 P 2 = Hqigi - 9202 ) , (40) 

where A = Tr\pipi—p 2 P 2 \- As pi / p 2 , we have 0 < A < 1 
according to Eq. Here, A > 0 holds as states of 
an optimal set are nonequal by definition. By means of 
linearity of the dynamical maps <i>t and homogeneity of 
the trace norm one finally obtains 

\\qi 0 i{t) - q 2 g 2 {t)\\i = j\\pipi{t) - P 2 P 2 {t)\\i (41) 

for all t > 0 where A“^ > 1. This shows that any in¬ 
crease of \\pipi{t) — P 2 P 2 {t)\\i is exceeded by the increase 
of \\qigi{t) - q 2 g 2 {t)\\i. Hence, {qi, g 1 , 2 } yields a non- 
Markovianity strictly larger than the set {pi,pip} con¬ 
tradicting its optimality. □ 

We define 5(77) to be the interior of the state space, i.e. 
the set of all quantum states p for which there is an e > 0 
such that all Hermitian operators X with unit trace sat¬ 
isfying Up — A||i < e belong to 5(77). Hence, all states 
in the interior have full rank. Note that 5(77) = 0 if 
dim77 = 00 which implies that the local representation of 
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the measure (37) introduced below is not available in in¬ 


finite dimensions. However, these quantum systems can 
frequently be accurately described by finite-dimensional 
Hilbert spaces. 

Based on the orthogonality of states of an optimal set 
we can also establish a local representation for the gener¬ 
alized measure as provided for the original definition in 
Ref. [2] . To this end, one first has to prove a character¬ 
ization of enclosing surfaces which are defined as follows: 
A set dU(p) C S{'H) not containing p G S{'H) is called 
an enclosing surface of p if and only if for any nonzero, 
Hermitian and traceless operator Y there exists a real 
number p > 0 such that 


p + 2pY G dU{p) . 


(42) 



Figure 1. Illustration of an enclosing surface dU (p) of an inner 
point p of the state space S{'H). 


Lemma 2. Let dU{p) he an enclosing surface and X a 
nonzero, Hermitian and indefinite operator. Then there 
exists a real number A > 0 with A|TrX| < 1 such that 

— {P+P - sgn(TrA)AA) G dU{p) , (43) 

P- 


where p± = ^(1 ± A|TrAi|) and 


sgn(a;) 


— 1 , if X < 0 

-fl , else 


(44) 


so that A = ^/(l-|-^|TrA|) > 0 holds and hence A|TrA| < 
1. For p± = 5(1 ± A|TrA|) one finally obtains 


1 _ X\TrX\ {®g’^(TrX) [(TrX)p - X] } 

= p-^sgn(TrX)X -|- — - —p 

P- P- 

= —{p+P - sgn{TrX)XX) , (48) 

P- 


which is Eq. (431. 


□ 


As proven in Appendix lemma can actually be 
augmented to show that the characterization (431 of an 
enclosing surface is equivalent to Eq. ^2 


Theorem 4. The generalized measure of quantum non- 
Markovianity admits a local representation, i.e. 


JV{^)= max / dt a{t,pi, Pi) (49) 

{Pi},P2&dU(pi) 7 ( t >0 


with 


d{t,pi,pi) = 


1 d 
IbiPi -P2P2II1 dt 


Ibl^t(Pl) -P2^t{p2) 


(50) 


where pi G S{'H) is any fixed inner point of the state 
space and dU{pi) refers to an arbitrary enclosing surface 
of Pi- 


We emphasize that a{t,pi,pi) is well-defined as 0 < 
IbiPi —P 2 P 2 II 1 for any state p 2 S dU{pi) and probability 
distribution {pi} as by definition of an enclosing surface 
we have p 2 b Pi (of- Eq. (§)• 


Proof. Let dU{pi) be an enclosing surface of pi. To 
prove that the corresponding local representation yields 
a value smaller than or equal to the original definition 
(391 one follows the lines of the proof of theorem Let 
P 2 G dU{pi) and a probability distribution {pi} be given. 
Now, according to lemmaj^there exist two state pi T p 2 
and a probability distribution {qi} such that 


PiPi - P2P2 = X{qiQi - P 2 P 2 ) , (51) 


Proof. Let dU (p) be an enclosing surface of p. Consider a 
nonzero, Hermitian, indefinite operator X. The operator 


where A = |biPi — P2P2II1 > 0. We then obtain 

^{t,Pi,Pi) = <j{t,qt,gi) , (52) 


Y = sgn(TrX) [(TrX)p - X] (45) 

defines a nonzero, Hermitian and traceless operator. By 
definition, there exists a real number p > 0 such that 

g= p+2pY edU{p) . (46) 


We then define a real number A by means of 

A 

^ “ 1 - A|TrX| 


(47) 


for all t > 0 due to linearity of <!)( and homogeneity of 
the trace norm and the derivative. We conclude that the 
right-hand side of Eq. (491 is smaller than or equal to 
Af{^) as defined in Eq. (39l. 

Conversely, let pi T p 2 be two orthogonal states 
and denote by {qi} a probability distribution. Then 
^ = Qigi ~ 92 P 2 defines a nonzero, Hermitian, indefi¬ 
nite operator. Thus, according to lemma there exists 
a real number A > 0 with A|TrA| < 1 such that 


p+pi - P-P 2 = cAA (53) 













where c = sgn(TrA) and p± = ^(1 ± A|TrA|) for some 
quantum state p 2 & dU{pi) of the enclosing surface. As 
|c| = 1 it follows that \\p+pi — P-P2II1 = A > 0. Linearity 
of the dynamical map and homogeneity of the trace norm 
and the derivative then yields 




(54) 


showing that the original definition (39) leads to a value 


which is smaller than or equal to the right-hand side of 
Eq. (491. This thus concludes the proof, that is, the max¬ 


imization over an enclosing surface with an information 
flux rescaled by the initial distinguishability reproduces 
the dynamics of the trace norm of orthogonal states. □ 


for 0 <t <ti with ti < T and 

+ ^-^[(^+p{t)(^--l{^-(^+.p{t)}] , (56) 

ii ti < t < T holds. Here, aj refers to the Pauli matrices 
and cr_|_(_) describes the usual raising (lowering) opera¬ 
tor with respect to the eigenstates |±) of (J 3 . Moreover, 
the rates satisfy j{t),b{t) > 0 for all t in their respec¬ 
tive domains so that the second phase of the process is 
neither CP- nor P-divisible which can be easily shown 
using theorem This dynamical process exhibits the 
essential feature of the generalized amplitude damping 
channel studied in Ref. [5T] . 


A careful reformulation of an enclosing surface thus al¬ 
lows to establish an equivalent local and universal repre¬ 
sentation of the generalized trace distance based measure 
as for the original definition m- This shows that also 
for the novel characterization, non-Markovianity is a uni¬ 
versal feature appearing everywhere in state space. Sim¬ 
ilarly, it suffices if Eq. (421 holds for exactly one p > 0 
as A is uniquely determined by this parameter given a 
nonzero, Hermitian and indefinite operator (cf. proof of 
lemma [^. In addition, also here, no assumption on the 
enclosing surface concerning for example the shape or 
the smoothness is needed implying a great benefit for the 
analytical, numerical and experimental determination of 
the generalized measure. 


B. Examples 


As already stated before the fundamental difference 
between the trace distance based measure for non- 
Markovianity and its generalization is that in the trace 
distance approach P-divisibility is only a sufficient cri¬ 
terion for Markovianity. This fact becomes particularly 
apparent when uniform translations of states are encoun¬ 
tered in the dynamics which do not describe positive 
maps but leave the trace distance unchanged |31j . How¬ 
ever, choosing unequal weights pi the Helstrom matrices 
are no longer invariant under such operations thus mak¬ 
ing it possible to detect the effect of such maps. 

To illustrate this property, we consider a two-level open 
quantum system, i.e. T-L = C^, undergoing a dynamics 
which, as depicted in Fig.[^ can be described in the Bloch 
sphere as an isotropic contraction followed by a transla¬ 
tion along the z-axis. The dynamics is analytically de¬ 
scribed by a time-local master equation (241 where the 
generator K-t obeys 


^tp(i) = ^ ^ Kp(t)crj - p{t)] , (55) 



t = 0 



t = T 



Figure 2. Bloch sphere representation of the action of the 
dynamical map 4>t corresponding to the time-local generator 
ICt defined by Eqs. (55 1 and (561 for t = 0,1 1 and T. 


Employing the Bloch vector representation for two- 
level systems the master equation is equivalently de¬ 
scribed by a differential equation for the Bloch vector. 


d _ f A{t)v{t) , 0 < t < ti < T 

dt 15(f) , ti <t <T 


(57) 


where A{t) = diag(— 7 (t), — 7 (t), — 7 (<)) and b{t) = 
(0, 0, 5(t))^. Hence, the process’ first phase corresponds 
to an isotropic contraction of the Bloch sphere Bi = 
{E| I El < 1} to the smaller sphere 5,. = {E| |E| < r} with 
radius 


r = exp 



e (0,1) , 


(58) 


which is clearly CP-divisible and, therefore, Markovian. 
On the other hand, the second phase describes a uni¬ 
form translation of the Bloch sphere along the z-axis, i.e. 
v{ti) I—>■ v{ti) + (0,0, a)"^ with 


a = / dtb{t) > 0 

J ti 


(59) 


where we thus have to require a < 1 — r in order to 
maintain positivity of the dynamical map (see Fig.[^. It 
is easily shown that this condition is also necessary and 
sufficient for complete positivity of the process. 
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The trace norm of the Helstrom matrix A = pipi — 
P 2 P 2 at time t for two quantum states pi ^2 = ^{'^2 + 'Vi ,2 ■ 
a), evolving according to the described dynamical map, 
is given by 

l|A(t)|ll = \ {\P1 -P2 + |w(i)| \ +\Pl-P2 - |w(t)| 1} , 

where w{t) = piVi(t) — p2V2{t). Clearly, if pi = P2 = 
1/2, then ||A(t)||i = ^\vi{t) — V 2 {t)\ showing that the 
unbiased case is unable to detect the non-Markovianity 
of the process resulting from the uniform translation. 
Without restriction we may assume pi > p 2 so that 


is valid for all 0 < t < and m ^ n which is con¬ 
dition (15l. Hence, the trace distance based measure 
and its generalization are equal to zero in this case while 
measures for non-Markovianity relying on CP-divisibility 
lainiisa are, of course, nonvanishing. This random uni¬ 
tary evolution [Ml [Ml illustrates the persisting and sig¬ 
nificant difference between the two major approaches for 
the characterization of quantum non-Markovianity. 


IV. CONCLUSIONS AND OUTLOOK 


l|A(t)||i 


P1-P2 , if Pi -P2> |w(t)| 
|u;(t)| , if Pi-P 2 < |w(t)| . 


(61) 


and, due to theorem [^ we may restrict ourselves to 
orthogonal states corresponding to antipodal unit vec¬ 
tors. That is, vi = —V 2 with Iwil = 1, which implies 

w{t) = Vi{t). 

If the probability distribution is such that pi — P 2 > 
1111(^)1 holds for all ti < t < T, then ||A(t)||i is a mono- 
tonically decreasing function of t. However, if pi — P 2 < 
|iTi(T)|, then 


[ dtcr(t) = |ui(T)| - (pi-P 2 ) > 0 , 
J a>0 


(62) 


indicating non-Markovianity. The change of the trace 
norm thus increases with decreasing difference pi — p 2 
which is bounded by pi — p 2 = r = |ui(<i)|. Calculating 
explicitly the change of the trace norm for pi — P 2 < r 
one finds 


||A(r)|li-||A(U)||i 

= + (pi - P 2 )a| - |ui(u)| , ( 63 ) 

which shows that the trace norm attains its maximal 
value 


||A(T)|U-|lA(ti)l|i=r-|a| , (64) 

if vi is parallel to a, i.e. vi = ca for some c > 0 , and 
Pi — P2 = T. Obviously, this is satisfied by the proba¬ 
bility distribution {pi} = {p± = i(l ± r)}. Unlike the 
original definition, the generalized trace distance based 
measure is thus able to capture non-Markovian dynam¬ 
ics arising from uniform translations contained in the dy¬ 
namical process. 

The dynamics generated by /Ct for 0 < t < U (cf. 
Eq.(l^) for any U > 0 provides also an example for 
a process which is not CP- but P-divisible. Choosing 
7i(t) = 72(t) = 1 and 73(t) = — tanh(<) as proposed 
in Ref. [M]? the dynamical process is not CP-divisible 
according to theorem In particular, there exists even 
no single interval for which CP-divisibility is restored as 
73(t) < 0 for all t > 0. However, the dynamics is always 
P-divisible since 
2 

|(m|tTi|n)p — tanh(<)|(m|CT 3 |n)p > 0 , (65) 

1=1 


We have introduced a generalization of the criterion 
of quantum non-Markovianity based on the flow of infor¬ 
mation. This novel characterization relies on the trace 
norm of Helstrom matrices which can also be interpreted 
as a measure for the distinguishability. By virtue of this 
property, the generalized measure still admits an inter¬ 
pretation as quantifier of an information backflow from 
the environment to the open quantum system. 

It is shown that the generalized criterion is equiva¬ 
lent to P-divisibility of the dynamical process which has 
an explicit connection to classical Markovian stochastic 
processes. That is, any rate equation obtained from a 
quantum master equation of a P-divisible process can be 
interpreted as Pauli master equation of a classical Markov 
process. However, the presented approach is more gen¬ 
eral since it can also be applied even when the notion of 
divisibility is ill-defined and can be experimentally tested. 

The experimental determination is substantially sim¬ 
plified by the derived mathematical representations of 
the generalized measure which are similar to those for 
the original definition. First, we have demonstrated that 
optimal initial states for non-Markovian dynamics must 
be orthogonal and, based on this result, we could finally 
establish a local representation for the measure. Hence, 
orthogonal states, corresponding to a maximal informa¬ 
tion content, exhibit maximal memory effects which can 
be revealed locally and anywhere in the quantum state 
space as provided by the local representation. 

An essential feature of the generalized approach to 
non-Markovianity in comparison with the original def¬ 
inition is its sensitivity to memory effects arising from 
uniform translations of states. To illustrate this, we 
constructed a dynamical process for a two-level sys¬ 
tem which comprises a uniform translation of the Bloch 
sphere. However, there exist dynamical processes which 
are not CP-divisible but P-divisible as we show in our 
second example which manifests the existing difference of 
the generalized definition with other approaches to non- 
Markovianity. 

We believe that the definition of non-Markovianity 
given in this paper is of great relevance for the study of 
memory effects in the field of complex quantum systems 
and quantum information due to the experimental acces¬ 
sibility and its clear-cut interpretation and connection to 
classical Markov processes. 
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Appendix A: Orthonormal bases and the maximally 
mixed state 


We prove the statement made in the main text (see Sec. 


IICI about the relation between the maximally mixed 


state and the eigenbases of the set of time-evolved states 
Imd)^ = {<i)t(/9)]p G S{TL)} for two-level systems, i.e. TL = 
C^. Note that Imd)^ defines a nonempty, convex and 
compact set for any finite-dimensional Hilbert space due 
to linearity of <i>t along with compactness and convexity 
of the state space. 


Lemma 3. Any orthonormal basis {['I'.i)} of PL defines 
the eigenbasis of a quantum state p(t) G Im^t if and only 
i/ 5 I 2 G Im<i)t. 


Proof. The ’if statement is clear from the fact that any 
orthonormal basis defines a resolution of identity. Hence, 
if ^l 2 G Imd)^ then any basis {Idfi)} defines at least 
the eigenbasis of the maximally mixed state proving the 
claim. It is clear that the same reasoning actually applies 
to any finite-dimensional Hilbert space. 

To show the reverse we employ the Hahn-Banach sep¬ 
aration theorem for normed vector spaces. Suppose that 
il 2 ^ Imd)^. Then there exists a real-valued, linear 
and continuous functional ipA separating the disjoint 
nonempty, convex and compact sets Imd)^ and {II 2 }. 
That is, we have 

'Pa{\'^ 2) <ini{q)A{p{t))\p{t) G Im$t} , (Al) 

where lpa{PC) = Tt{AX} for some Hermitian opera¬ 
tor A due to Riesz representation theorem. By means 
of the transformation A 1 —>■ A — TrA one may also as¬ 
sume that the operator A is traceless. The set {X\X = 
Al, Tr{A} = 1, ipa{X) = 0} thus describes a hyperplane 
which intersects the maximally mixed state and separates 
it from Im<i>t. 
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Figure 3. Two-dimensional cut of a hyperplane in the Bloch 
representation separating the two disjoint nonempty, convex 
and compact sets {II 2 } and Im<l>t. As the plane intersects the 
maximally mixed state its intersection with the state space 
comprise an uncountable set of pure and orthogonal states 
defining orthonormal bases of C^. 


Employing the Bloch representation this hyperplane 
defines an ordinary plane in which contains the origin 
and therefore intersects the surface of the Bloch sphere 
in a unit circle. As pure, orthogonal states on cor¬ 
respond to antipodal points on the surface of the Bloch 
sphere we thus conclude that the (hyper)plane contains 
an uncountable set of orthonormal bases which do not de¬ 
fine eigenbases of quantum states p{t) G Im^t (cf. Fig. 

i- 

□ 


Appendix B: Characterizations of enclosing surfaces 

In this section we show that the relation for enclosing 
surfaces in terms of nonzero, Hermitian and indefinite 
operators derived in lemmaj^is equivalent to the original 
definition, i.e. 

Lemma 4. A set dU (p) defines an enclosing surface of 
p if and only if for any nonzero, Hermitian, indefinite 
operator X there exist a real number A > 0 with A|TrA| < 


1 such that 


— {P+P - sgn{TTX)XX) e dU{p) , 
P- 

where p± = |(1 ± A|TrA|) and 


sgn(a;) 


— 1 , if X < 0 

-1-1 , else 


(Bl) 


(B2) 


Proof. That any enclosing surface obeys a characteri¬ 
zation in terms of nonzero, Hermitian and indefinite 
operators has already been proven in lemma Con¬ 
verse^ suppose the states in dU(p) are characterized by 
Eq. ( |B1[ ). Hence, for a nonzero, indefinite and Hermi¬ 
tian operator X there exists a real number A > 0 with 
A|TrA| < 1 such that 

p = — {p+p — sgn(TrA)AX) G dU{p) (B3) 

Pt 

where p± = ^(1 ± A|TrA|). Now, consider the map 

0p(A) = sgn(TrA) [(TrA)p - A] (B4) 

defined on the set of nonzero, Hermitian and indefi¬ 
nite operators. The operator Y = Qp{X) represents a 
traceless, Hermitian operator and E = 0 if and only if 
(TrA)p = X contradicting that X is indefinite. Thus, 
we have Y ^ 0 and one finds 

P = ~ {p+p - - Y]} 

p+-A|TrA| ^^ A^, 

P- P- 

= p + 2pY (B5) 

where p. = A/(2p_) > 0. It remains to show that 0p 
defines a surjection on the set of nonzero, Hermitian and 
traceless operators. This is obviously true as any trace¬ 
less, nonzero, Hermitian operator Y is necessarily indef¬ 
inite and we have 0p(E) = Y. Hence, the set dU{p) is 
indeed an enclosing surface. □ 



